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Abstract 

We demonstrate two periodic or quasi-periodic generalizations of the Chaplygin gas (CG) type 
models to explain the origins of dark energy as well as dark matter by using the Weierstrass p{t), 



^SJ ' a{t) and C,{t) functions with two periods being infinite. If the universe can evolve periodically, a 

^O ' non-singular universe can be realized. Furthermore, we examine the cosmological evolution and 

^ , 

Cn . nature of the equation of state (EoS) of dark energy in the Friedmann-Lemaitre-Robertson- Walker 

'NT ■ 

ff^ \ cosmology. It is explicitly illustrated that there exist three type models in which the universe always 

o 

CN ' stays in the non-phantom (quintessence) phase, whereas it always evolves in the phantom phase, 

or the crossing of the phantom divide can be realized. The scalar fields and the corresponding 
/\ , potentials are also analyzed for different types of models. 
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I. INTRODUCTION 



Inflation in the early universe has been confirmed by the recent observations of cosmic 



M 



microwave background (CMB) radiation [11-141 • ■'■^ addition, the accelerated expansion of the 
current universe has also been suggested by recent observations, e.g.. Type la Supernovae 5|, 
6|, CMB radiation \MiM^ the large scale structure LSS JTI, |8|, baryon acoustic oscillations 
(BAO) |9|, and weak lensing 10|. To explain such a cosmic acceleration, one provides the 
existence of so-called dark energy in the framework of general relativity (for reviews, see. 



e.g., 
e.g., 



and 



11-18[), or one supposes that gravity is modified on the large scale (for reviews, see, 

id-M). 

In the expansion history of the universe, there exist two singularities. One is a Big Bang 
singularity. The other is the finite-time future singularities 27H74|. which occurs at the last 
stage of the universe filled with dark energy, or a Big Crunch singularity. To avoid these 
singularities, various cosmological scenarios have been proposed, e.g., the cychc universe 
85| (for a reference in a d iffere nt context, see 86|), the ekpyrotic scenario 
the bouncing universe 90Nl0l| . Furthermore, related to t he cyclic universe, the (trefoil and 
figure-eight) knot universe has been investigated in Ref. 102Nl05| . In addition, motivated 
by the studies on the role of applying the Weierstrass p{t), C (t) and a (t)-functions and 
the Jacobian elliptic functions to astrophysics and cosmology 106l - lll0| . the equation of 
state (EoS) for the cyclic universes in the homogeneous and isotropic Friedmann-Lemaitre- 
Robertson- Walker (FLRW) spacetime has been reconstructed by using the Weierstrass and 
Jacobian elliptic functions in Ref. 111 |. 
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71j . with the 



In this paper, based on the reconstruction method in Refs. [19|, 
Weierstrass p(t)-function, we examine the cosmological evolution of the EoS for dark energy 
in FLRW cosmology. In particular, it is shown that two periodic generalized Chaplygin gas 
(GCG) type models for dark energy can be reconstructed. To accou nt for t he o rigins of 
dark energy as well as dark matter with a fiuid, the original CG 112| . GCG 113| and the 
modified CG (MGC) 114J . Ill6| have been explored. We mention that the reconstruction 
of periodic cosmologies have widely been stud ied. Especially, the reconstruction of periodic 
EoS has been investigated, e.g., in Ref. 115| . In this reference, with an inhomogeneous 
EoS for dark energy fiuid, it has been demonstrated that an oscillating universe can occur. 
Also, the Hubble parameter with a periodic behavior can realize both infiation in the early 



universe and the late-time cosmic acceleration under the same mechanism in a unified man- 
ner. In addition, it has been verified that a coupling between dark energy fiuid, which has a 
homogeneous and constant EoS, and matter, can present a periodic behavior of the universe. 
Furthermore, as several theoretical issues in the universe with its oscillatory behavior, the 
phantom phase and finite-time future singularities have been investigated. A scalar-tensor 
description of the oscillating universe has also been explored. As stated, there exist various 
theoretical subjects in the periodic cosmological evolution of the universe. The essential 
property of the Weierstrass functions is to have two periods mi and ?7i2- Hence, the period- 
ical and quasiperiodical models which we show in Sec. 3 are periodical or quasiperiodical in 
terms of the energy density p. In addition, these models with the periods rrii and 7712 being 
infinite are reduced to the Chaplygin gas models (as is seen from the formulae in Eq. fl3.ll) ). 
Thus, the reconstruction procedure of these models corresponds to two periodic or quasi pe- 
riodical generalizations of the CG models. This justifies the use of the Weierstrass functions 
in cosmological models. Furthermore, the models given in Sec. 4 are periodic on the cosmic 
time t (which is a dimensionless quantity in our analysis), although the models in Sec. 3 
are periodic/quasiperiodic on p. The periodicity of the cosmological evolution comes from 
the periodic nature of the Weierstrass functions. Also, the periodic/quasiperiodic models in 
Sec. 3 are singular at p = 0, whereas those in Sec. 4 are singular at t = 0. If the periodic 
evolution of the universe can be realized, various scenarios to avoid cosmological singularities 
can be constructed. This is the important cosmological motivation to obtain such periodical 
solutions. Moreover, we explicitly demonstrated that there exist three type models in which 
(i) the universe always stays in the non-phantom (quintessence) phase, (ii) it always evolves 
in the phantom phase, and (iii) the crossing of the phantom divide can be realized. I t ha s 
recently been shown that these three cases have also been realized in non-local gravity 117 |. 



It is also interesting to remark that according to the analysis of recent c osmologi cal obser- 



vational data, in the past the crossing of the phantom divide occurred [I18l - ll22l ]. We use 
the units of the gravitational constant SirG = c = 1 with G and c being the gravitational 
constant and the seed of light. 

The paper is organized as follows. In Sec. 2, we show the basic equations in the FLRW 
background and briefiy give the Chaplygin gas type models. In Sec. 3, we study periodical 
and quasi-periodical GCG type models. In Sec. 4, we demonstrate other two periodical 
FLRW models. Finally, several conclusions are presented in Sec. 5. 



II. BRIEF REVIEW OF THE CG TYPE MODELS 



In this section, we briefly explain the significant features of the CG type models for the 
spatially flat homogeneous and isotropic FLRW universe. The action describing general 
relativity and matter is given by S* = J y/—gd^x{R + Lm)-, where R is the scalar curvature 
and Lm is the matter Lagrangian. We take the flat FLRW spacetime with the metric, 
ds^ = —dt^ + a^{t) {dr"^ + r^rffi^). Here, a{t) is the scale factor and dQ"^ is the metric of 
2-dimensional sphere with unit radius. We note that in this paper, time (t) is considered 
to be a dimensionless quantity. In the flat FLRW background, from the above action we 
obtain the gravitational field equations 

(2.1) 

(2.2) 

Here, the Hubble parameter is defined by if = d/a and a dot denotes the time derivative of 
d/dt. By using these equations, we have the expressions of the energy density p = 3 (d/a) 
and pressure p = —2 (d/a) — (d/a) . 

Next, we explore the CG type models 
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116| . The GCG model has been con- 



structed in order to account for both the origins of dark energy and dar k matter with using 
a single fluid. The equation of state (EoS) of the GCG is given by 113 1 

P = -^^ (2.3) 

where Ci(> 0) is a posit ive c onstant and a is a constant. If we take a = 1, Eq. (12. 3 p describes 
the original CG model 112|. From Eq. (12. 3p and the continuity equation p = —3H{p + p), 
we obtain 



P 



Ci + 



Co 



l/{l+a) 



(2.4) 



a3(i+")_ 

where C2 is a constant of integration. From Eq. (12. 4p . we find the asymptotic behaviors 
of p that in the early universe a ^ 1, p ~ C2a~^, whereas in the late universe a ^ 1, 
p ~ Ci " . Thus, in the early universe, the energy density behaves as p oc a~^, which 
is the same as non-relativistic matter such as dark matter. On the other hand, in the late 
universe the energy density becomes a constant as p —^ C-^^ ° . This means that it can 
play a role of dark energy. As a consequence, the GCG model can explain the origin of dark 
energy as well as dark matter simultaneously. 



In addition, the MCG has been proposed in Refs. 114l . Ill6| . The EoS is given by 



P = Cap - ^ , (2.5) 



where C3 and C4(> 0) are constants. 



III. PERIODICAL AND QUASI-PERIODICAL GCG TYPE MODELS 

In this section, we examine the periodical and quasi-periodical GCG type models by using 
the Weierstrass functions, the so-called MG - i models^. Properties of elliptic functions 
inform us that the MG - XXV, MG - XXVI and MG - XXIV models (as the MG - XXIII 
and MG - XXVI models) are some generalizations of the CG type models due to the following 
degenerate cases of som e elli ptic and related functions as rrii = 00 and 7712 = 00, where rrii 



and 7712 are two periods 123| : 



cr(x) = X , ({x) = X , p{x) = X , a7n{x) = X . (3.1) 

Here, (72 = fi's = and ei = 62 = 63 = 0, where g2 and g^ are the Weierstrass invariants. 

The physical motivation to examine the series of the MG - i gas is as follows. These 
models can realize the cosmological evolution of the GCG type models with the periodical 
and quasi-periodical behaviors, which depends on the models. These models are expressed 
with the Weierstrass functions and hence various behaviors of the cosmic expansion history 
with periodicity and/or quasi-periodicity can be realized. Thus, these models can present 
novel cosmological scenarios without a Big Bang singularity in the early universe and the 
finite-time future singularities or a Big Crunch singularity, such as the cyclic universe, the 
ekpyrotic scenario and the bouncing universe. 



1 Note that the meaning of the so-called Myrzakulov Gas (MG) - i {i = XII, XIII, XIV, XV, XVI, XVII, 
XVIII, XIX, XX, XXXIII, XXI XX II, XXII, XXIII, XXVII) is the model of some gases/fluid, which is 



the notation used in Refs. 



102 



105| 



A. Periodical generalizations 



1. MG- XXI model 



One of the most int eresting examples of gases is the MG - XXI model which has the 



following EoS [102 



-1051 



V 



-B[p{p)r , (3.2) 

where B{> 0) is a positive constant. By using the de gene rate case of the function p(p) 
in (13. ip . we can show that the well-known CG model 112] p = —B/p, which is equal to 
Eq. f l2.4p with B = Ci and a = 1, is particular case of the MG - XXI model in Eq. f l3.2p . 
The parameter of the EoS u for our model is given by 



P 

oj = — 

P 



P 



(3.3) 



2. MG - XXII model 



One of the two periodical generalizations of the GCG is given by 102l - ll05 | 



p = -B[p{p)] 



0.5a 



(3.4) 



In fact, its degenerate case is the GCG |113J p = —B/p°', which is equivalent to Eq. (12. 3p 
with B = Ci. For this model, the parameter of the EoS looks like 



u 



(3.5) 



3. MG - XXIII model 



!^ext, we present one of the two periodical generalizations of the MCG. Its EoS reads 102 - 

105| 

p = Ap-S[p(p)^^ (3.6) 



where A is a constant. The corresponding parameter of the EoS is 

B[pip)r^ 



uj = A 



P 



(3.7) 



4. MG - XXIV model 

We now give a more general form of two periodical generalizations of the MCG. Its EoS 
is described as 

P = A[pipr'-' - B[p{p)r-'- . (3.8) 

The parameter of the EoS for the model is written by 

,, _ A[p{p)r-' i?b(p)]°-^" .gg. 

P P ' 

Again, by using the degenerate properties of the elliptic functions, we can demonstrate that 
this model is reduced to MCG. 

B. Quasi-periodical generalizations 

In the preceding subsection, we have considered two periodical generalizations of CO type 
models. In this subsection, we study quasi-periodical models. 

1. MG- XXV model 

One of the quasi-periodical models, the so-called MG - XXV model, is given by 

where A and B are constants and a{p) is the Weierstrass a-function. As the cr-function 
degenerates according to equations in (13.11) . in this case Eq. (I3.10p becomes the MCG in 
Eq. (12. 5p . The corresponding parameter of the EoS is expressed as 

Aaip) B 

u = -— — — . (3.11) 

2. MG - XXVI model 

Our next quasi-periodical model is given by 

P = ^-5[C(P)]", (3.12) 



where ({x) is the Weierstrass C(x)-function. It is the MG - XXVI model. For the degenerate 
case in Eq. (13. ip . this model is also reduced to the MCG. Its parameter of the EoS takes 
the form 

c.--^-^^^. (3.13) 

pC(p) p 

3. MG - XXVII model 

We explore the MG - XXVII model. For this model, the EoS reads 

p = Aam{p) - 5[am(p)]-" , (3.14) 

where am{x) is the Jacobi amplitude {am{x)) function and a is a constant. In case of the 
degeneration in Eq. (13. ip . this model recovers the MCG. The parameter of the EoS is written 

by 

Aamjp) g[am(p)]~" 

uj = , (3.15) 

P P 

It is significant to emphasize that (a) if we substitute a{t) in Eq. (13. ip into Eq. (I3.10p . (b) 
if we use ({t) in Eq. (O) and Eq. fl312D . (c) if we combine p{t) in Eq. (O) with Eq. fl314l) . 
then we obtain the MGC in Eq. (12. 5p . As a result, in the limit of rrii = oo and m^ = cxp , 



114 



116 | 



the MG - XXV, MG - XXVI and MG - XXIV models are reduced to the MCG 
This point is the most important and novel observataion in this work. 

In the limit of the small energy density p — )■ as well as p — ?■ oo, the behaviors of the EoS 
for the universe in the MG - XXI, MG - XXIV, MG - XXV, MG - XXVI and MG - XXVII 
models asymptotically approach those in the CG model. On the other hand, in the middle 
regime of p, since the EoS for the universe in the MG - XXI, MG - XXIV, MG - XXV, MG 
- XXVI and MG - XXVII models is described by using elliptic functions with a periodic or 
quasi-periodic property, the EoS for the universe expresses also periodic or quasi-periodic 
behaviors. 

From the above considerations, the cosmological evolution of the universe is described 
as follows. First, the energy responsible for inflation would be released to radiation (i.e., 
relativistic matter) through a reheating process and the universe enter the radiation domi- 
nated stage. Here, the concrete mechanism for both inflation and the reheating stage is not 
specified. After that, as the universe expands, its temperature decreases with proportional 
to a~^, and the matter (i.e., non- relativistic matter) dominated stage appears. This can 

8 



be seen in our models in the limit of p — )■ oo, namely uj asymptotically approaches zero, 
which corresponds to the EoS of the dust. Finally, the universe becomes the dark energy 
dominated stage. This can also be understood in the limit of p — )■ 0, where u < —1/3. 
Thus, it is considered that the cosmological evolution of the universe can be realized in our 
models. 

IV. OTHER TWO PERIODICAL FLRW MODELS 

The EoS for dark energy is one of the most significant cosmological quantities. In this 
paper, we concentrate on the evolution of the EoS for dark energy. In the FLRW spacetime. 



the effective EoS for the universe is given by [l9|, l20| Weg = Pcs/PcS = — 1 — 2H/ (SH^). 
Here, pcs and pcs can be considered to the total energy density and pressure of the universe, 
respectively. Since we examine the dark energy dominated stage, the energy density pde 
and pressure pde of dark energy can be regarded as pde ~ PcS and pde ~ PcR- As a result, 
we find wde ~ <^e{i- 

In addition, we represent pde and Pde as p in Eq. (12. 2p and p in Eq. (12. ip . respectively. 
In the non-phantom (quintessence) phase, H < and hence Ucs > — 1 , which is the non- 
phantom (quintessence) phase, while in the phantom phase, H > and therefore Wcflf < — 1. 
li H = 0, Ucs = —1, which is the case that dark energy is the cosmological constant. 

As a qualitative criterion to constrain the models, we examine the evolution of the EoS u 
of a fluid corresponding to dark energy. If uj is always less than —1, the universe stays in the 
phantom phase in all the cosmic evolution history. This case is clearly inconsistent with the 
standard cosmological evolution and hence it can be ruled out. On the other hand, if u is 
always larger than —1 or it crosses the line of —1, these cases are not ruled out, namely these 
models may have the possibility to realize the standard evolution history of the universe. 

In this section, we study new two periodical FLRW models. These models are expressed 
by using the Weierstrass p(t)-function, which as well known satisfies the following ordinary 



differential equations [123 1 



p2(t)= V(t)_^2p(t)-^3, 

f3(t) = 6p2(t)- 0.5^2, 
P{t) = 12p{t)p{t), 
'P{t) = 120p\t)-18g2pit)-12g3 



(4.1) 
(4.2) 
(4.3) 
(4.4) 



where p( t 



Refs. 



19 



20 



dp(t) Idt and so on. In what follows, by using the reconstruction method in 
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7l|, and the Weierstrass p(t)-function, for ten models [the MG 



models {i = XII, XIII, XIV, XV, XVI, XVII, XVIII, XIX, XX, XXXIII)], we reconstruct 
the EoS for dark energy and explore its cosmological evolution in FLRW cosmology. 

A. MG XII model 



We suppose that the Hubble parameter is represented as 



H = p{t) . 



From this expression, the scale factor becomes 



(4.5) 



a{t) = aoe"^(*) 



(4.6) 



where ao(> 0) is a positive constant and ((t) is the Weierstrass C(^)-function. Then, Eq. ( 14. 5p 
and the gravitational field equations (12. ip and (12. 2p lead to the parametric EoS 



p = -2^4p^{t)-g2p{t)~g3 - 3p\t) , 
p = 3p\t). 

By using Eqs. (14. 7p and (14.80 . we see that the EoS parameter is written as 



(4.7) 
(4.8) 



^P^,_ 2^/4p%t) - g^pjt) - gs 
^ p 3p\t) 



(4.9) 



In order to describe our models in terms of the scalar field theory, we introduce a scalar 
field and its self-interacting potential V{(j)). The Lagrangian for the scalar field theory is 
given by (see, e.g., |l9l. |20|) 



10 



Thus, this scalar is corresponding to a phantom one with u < —1, which can be seen in 
Eq. f l4.9p . The energy momentum tensor of the scalar field (f){t) is identical to a fiuid with 
the energy density p^ and pressure p^ given by 

p<^ = -0.502 + 1^(0) =p, ^4^^^) 

p^ = -0.5<P^-V{cj))=p. (4.12) 

By using these expressions, we find 

-0' = P + P, (4.13) 

V{^) = 0.5{p-p). (4.14) 

In addition, it follows from Eqs. f l4.13p and fl4.14p that the scalar field </> and self- 
interacting potential V{(j)) are written as 

<P = t^j ^J^Ap\t) - g2p{t)^3dt , (4.15) 

V = 3p2 (t) + ^4p3(^t)-g2p{t)-gs. (4.16) 

In Fig. m we show the cosmological evolution of EoS w as a function of t for p (t, 1, 1), 
i.e., the model parameters of the Weierstrass invariants oi g2 = 1 and g^ = 1. From Fig. [H 
we see that the universe always stays in the phantom phase {u < —1). Hence, this model is 
ruled out. Furthermore, we find the two periodic oscillatory behavior of u. 

We remark a point in terms of the expression of V. In the above procedure, first the form 
of the scale factor a = a{t) or the Hubble parameter H = H{t) is supposed. Next, from 
this form we obtain the pressure p = p{t) and the energy density p = p{t). On the other 
hand, in the description of the scalar field theory, the scalar field = 0(t) and its potential 
V = V{(j)) are expressed with p = p(t) and p = p{t). Accordingly, V = V{(j)) = V{(f){t)). 
This means that V is written as a function of cosmic time t. Thus, in principle, if a, H, p 
and p are inversely solved in terms of t, V can also be described by the expressions of a, H, 
p and p. 

B. MG XIII model 

We express the Hubble parameter as 

H = p{t) . (4.17) 

11 




FIG. 1: The EoS uj in Eq. (I4.9P as a function of t for p (t, 1, 1), i.e., the model parameters of the 
Weierstrass invariants of (72 = 1 and (73 = 1. The line of w = — 1 is also plotted. 



From this expression, the scale factor is given by 



a{t) = aoe 



Pit) 



(4.18) 



By combining Eq. (I4.17P with the gravitational field equations (12. ip and (12.21) . we find that 
the parametric EoS is written as 



p = -12p^{t) - 12p\t) + 3g2p{t) +92 + 3(73 
p = 3 {4p%t) - g2p{t) - gs) . 



(4.19) 
(4.20) 



It follows from Eqs. (14.191) and (I4.20p that the EoS parameter is given by 

12p3(t) + I2p\t) - 3g2p{t) -g2- 3gs 



u 



(4.21) 



3{Ap^{t)-g2p{t)-gs) 

By using Eqs. (I4.13P and (I4.14p . we obtain the expressions of the scalar field (p and self- 
interacting potential V{(f>) 



-I2p^ it) + g2dt , 
V = I2p' (t) + 6p2 it) - 3g2P (t) - 3(73 - 0.5^2 • 



(4.22) 
(4.23) 



In Fig. 121 we depict the cosmological evolution of EoS w as a function of t for p [t, 1, 1) 
From Fig. |2l we see that the universe always evolves within the phantom phase {u < — 1) 
Consequently, this model is ruled out. 
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FIG. 2: The EoS w in Eq. (I4.2ip as a function of t. Legend is the same as Fig. [H 



C. MG XIV model 



We provide that the Hubble parameter is described as 



H = Pit) 



so that the scale factor can read as 



a{t) = aoe 



P(i) 



(4.24) 



(4.25) 



From Eq. (I4.24p and the gravitational field equations (12. ip and (12. 2p . we have the parametric 
EoS 



p = -24p(t) v/4p3(t) - gMt) -93-5 {Qp'{t) - O.^g^) , 
p = 3{6p\t) -0.592)' . 

Hence, the EoS parameter is written as 

, 8p{t)^4p^{t)-g2p{t)-gs 

U! = —1 o • 

{6p^{t) - 0.5^2) 
With Eqs. f liTTSj) and f04|) . we have 

(j) = 2l^ j ^/p (t) v/4p3 (t) _ g^p (t) - g,dt , 

V (0) = 3 (6p2 (t) _ 0.5^72)' + I2p (t) ^4p^t)-g2p{t)-g, 



(4.26) 
(4.27) 

(4.28) 

(4.29) 
(4.30) 



In Fig. [31 we illustrate the cosmological evolution of EoS w as a function of t for p(t, 1,1). 
From Fig. [3l we see that the universe always stays in the non-phantom (quintessence) phase 

(a;>-l). 
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FIG. 3: The EoS w in Eq. (|4.28p as a function of t. Legend is the same as Fig. [TJ 



D. MG XV model 



We take the form of the Hubble parameter as 

H = P{t) . 
From this form, the corresponding scale factor is given by 

a{t) = aoe^(*) . 



(4.31) 



(4.32) 



With Eq. (I4.3ip and the gravitational field equations (I4.13P and (I4.14p . we obtain the para- 
metric EoS 



P 



-12 {20p'{t) - 3g2p{t) - 2g,) - 432p\t) {4p'{t) - g^pit) - g-,) , (4.33) 



p = A32p{tf {Ap\t) - g2p{t) - g:,) . 
It follows from Eqs. (I4.33P and (I4.34p that the parametric EoS is given by 



u 



2Qp\t) - 3g2p{t) - 2^3 



(4.34) 



(4.35) 



(4.36) 
y (0) = 6 {20p%t) - 3g2p{t) - 2g,) + 432^^ (t) [Ap^ (t) - g^p (t) - g,) . (4.37) 

In Fig. m we display the cosmological evolution of EoS w as a function of t for p (t, 1, 1). 
From Fig. HI we understand that the universe always evolves within the phantom phase 
{uj < —1). Thus, this model is ruled out. 



36pity (Ap^it) - gMt) - 93) 
By using the formulae in Eqs. (14.131) and (14.140 . we acquire 

</) = 2z V3 [ v/20p3 (t) _ 3g2P (t) - 2g3dt , 
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FIG. 4: The EoS w in Eq. (I4.35P as a function of t. Legend is the same as Fig. [H 
E. MG XVI model 

We assume that the Hubble parameter is written as 

H = "p{t). 
In this case, the scale factor becomes 

a{t) = aoe'^'W . 
Using Eq. fl4.38p and the gravitational field equations (12.11) and (12. 2p yields 

p = -36 (20p2(t) _ g,) ^4p3(^t) - g2p{t) - 93 - 3 {l20p^{t) - ISg^pit) - Ug^)' , 
P = 3{l20p'{t)-18g2p{t)-12g3)\ 

Hence, the parametric EoS is expressed as 

. 12 {20p\t) - g^) V4p3(t) - g^pit) - g-s 
(120p3(t) - 18(72p(t) - 12(73)' 
From Eqs. (I4.13P and (I4.14p . we acquire 

= 6^ y V'(20p2(t) - ^2) v/4p3 it) - g^p (tf^^dt , 

V^ = 18 {20p\t) - (72) V'^PHt)-92p{t)-g2 + 3 (I20p3 (t) - I8g,p (t) - Ug,)' . 



(4.38) 

(4.39) 

(4.40) 
(4.41) 

(4.42) 

(4.43) 
(4.44) 



In Fig. we plot the cosmological evolution of EoS w as a function of t for p (t, 1, 1). 
From Fig. 0, we find that the universe always evolves within the phantom phase {u < —1). 
Therefore, this model is ruled out. In addition, we clearly see the two periodic oscillating 
evolution of oj. 
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FIG. 5: The EoS w in Eq. (|4.42p as a function of t. Legend is the same as Fig. [TJ 



F. MG XVII model 



We take the scale factor as 

« (t) = p{t) ■ 
From this expression, the Hubble parameter becomes 



H 



v/4p3 (t) - g2p (t) - gs 
Pit) 



(4.45) 



(4.46) 



Equation (I4.46P and the gravitational field equations (12. ip and (12. 2 p present the parametric 

EoS 

16pHt)-2g2p{t)~g3 



P 



p = 3 



pHt) 

4p' (t) - g2p (t) - 3^3 
PHt) 



Furthermore, the EoS parameter is represented as 

16p^t)-2g2p{t)-gs 



u 



12p3(t)-3^2p(t)-3^3 
From Eqs. (I4.13P and ( ]4.14p . we have 

v/V(t) + ^2p(t) + 2(73 



V 



Pit) 
28p2 (t) - 5g2P it) - 4(73 



dt , 



(4.47) 
(4.48) 

(4.49) 

(4.50) 
(4.51) 



2p2 (t) 

In Fig. El we plot the cosmological evolution of EoS cj as a function oft for p (t, 1, 1). From 
Fig. [6], we see that the universe always stays in the phantom phase (u; < —1). Accordingly, 
this model is ruled out. 
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FIG. 6: The EoS w in Eq. (I4.49P as a function of t. Legend is the same as Fig. [H 
G. MG XVIII model 

We describe the scale factor as 

a = pit) . 

It follows from this description that the Hubble parameter is given by 

6p2 (t) - 0.5^2 



H 



V4p3 it) - g2p (t) - 93 



H. MG XIX model 



We consider that the scale factor is given by 

a = p{t) , 
so that the Hubble parameter can read 



H 



6p2 (t) _ 0.5^2 



I. MG XX model 



We suppose the following form of the scale factor 

a= P{t). 

From this representation, the Hubble parameter is written by 

20p3(t)_3^2p(t)-2(73 



H 



Ip it) ^4p^t)-g2p{t)-g3 



(4.52) 



(4.53) 



(4.54) 



(4.55) 



(4.56) 



(4.57) 
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J. MG XXXIII model 

We suppose that the scale factor takes the following form 

a = "p{t). (4.58) 

It follows from this representation that the Hubble parameter is described by 

jj ^ S{g2-20p'{t)){g, + g,p^t)-Ap^t)) 

2g3 + 3g2p{t)-20p^t) ' ^ ' ^ 

It is important to mention that a quintom model (where there must exist both canonical 
and non-canonical scalar fields in order that the crossing of the phantom divide can occur) 
with single canonical scalar fi eld cannot be reconstructed, in which the crossing of the 



phantom divide cannot occur |l24j . Since in MG - XVIII, XIX, XX and XXXIII models 
explained in Sections 4.7, 4.8, 4.9 and 4.10, respectively, if we reconstruct the corresponding 
single scalar field model, the crossing of the phantom divide can be realized. Thus, this 
implies that for these models, such a reconstruction is not physical but just mathematical 
procedure. 

The physical interpretation of our results is considered as follows. By reconstructing the 
expansion behavior of the universe, i.e., the Hubble parameter, with the Weierstrass p{t)- 
function, which has a periodic property, the periodic behavior of the EoS for the universe can 
be realized. Such a procedure can be applied to scenarios to avoid cosmological singularities 
and eventually a non-singular universe can be realized. 

It is also remarkable to note that in terms of all the numerical calculations in Figs. [1]- 
El the qualitative results do not strongly depend on the initial conditions and the model 
parameters such as the Weierstrass p(t)-function. 

Moreover, for all the models, the scale factor a has to be positive and the Hubble pa- 
rameter H should be real. From Eqs. (gSD, (1418|) . ( K25\i . (1432]) and ( lOOj) . a is always 
positive. Also, it follows from Eqs. (I4.45p . (I4.52p . (I4.54p . (I4.56p and (I4.58p that a is written 
by p{t) or its time derivatives. By using the Laurent expansion of p{t) around t = as 
Pit) = l/t' + ig2/20)t' + igs/28)t^+Oit^),weRndpit) = -2/t^+ig2/10)t+igs/7)t^+O{t'), 
pit) = e/t^ + (72/10 + (3(73/7) t^ + 0{t^), Pit) = -24/t5 + (6(73/7) t + Oit'), and p" (t) = 
120/t^ + 6^3/7 + Oit^). Thus, for the cases of a = pit) in Eq. (|432|) and a = Pit) in 
Eq. (I4.56p . around t = 0, the expression of a is not well defined because the value of a would 
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become negative around t = 0. In this sense, around t = 0, the MG - XVIII and MG - XX 
models cannot be used. On the other hand, since the argument of the Weierstrass functions 
used in this paper is the cosmic time t, which is real, all the values of the Weierstrass func- 
tions are real. In all the models, H is described by the Weierstrass functions and hence H 
would be real. 

Furthermore, we state that if the evolution of u is displayed as a function of the redshift 
z = flp/a— 1 with ttp the present value of a, the direction of the cosmic time t to go by becomes 
opposite to that for u to be shown as a function of t in Figs. [THSl It would be considered that 
the main qualitative difference is only this point and hence other cosmological consequences 
do not change. 

In addition, it is interesting to note that the models examined in this work may solve, not 
only the dark energy paradigm, but also the horizon and flatness problems. In other words, 
these models may produce a kind of inflationary epoch or perhaps a contraction phase as 
in the ekpyrotic scenario. In order for inflation or the ekpyrotic scenario to be realized 
realistically, it is necessary to consider the way of connect the inflationary stage and the 
dark energy dominated stage, i.e., the realization of the reheating stage. This is a crucial 
point for these models to be realistic inflation models or the ekpyrotic scenario. 

Moreover, it should be cautioned that in Figs. ISHU apparently it looks the EoS u diverges. 
However, the reason is just a way of plotting and therefore there is no divergence in terms 
of the EoS in all of the figures. 

Finally, we mention the issue of existence of a ghost, namely, the instability for the case 
of the crossing of the phantom divide. For a scalar field theory, it is known that if the 
crossing of the phantom divide happens in the FLRW universe, there would appear a ghost. 
In Sec. IV, we have presented the interpretations of our models in the framework of a scalar 
field theory. However, there still exists the possibility that our models could be regarded as 
a more complicat ed theory , e.g., which is described b y the non -linear kinetic terms such as 



fc-essence models 
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5Hl30[ and the Galileon models 131nl35|. These investigations might 
be one of the interesting future works in our approach to the periodic and quasi-periodic 
generalizations of the Chaplygin gas type models. 
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V. CONCLUSIONS 

In the present paper, we have reconstructed periodic and quasi-periodic generahzations 
of the Chaplygin gas type models by using the Weierstrass pit), a{t) and C(^)-elhptic func- 
tions. We have explored the cosmological evolution of the EoS for dark energy in FLRW 
spacetime. In particular, we have shown that by using the degenerate properties of the 
elliptic functions, the MCG models can be recovered. This is one of the most important 
cosmological ingredients obtained in this paper. This result implies that by plugging the 
reconstruction method of the expansion history of the universe with the Weierstrass elliptic 
functions, we can acquire the MCG models, which have a potential to reveal the properties 
of both dark energy and dark matter. In other words, the procedure demonstrated in this 
paper can lead to a preferable cosmological model by starting with the mathematical special 
functions. 

It is meaningful to summarize the following significant points, (a) The Weierstrass func- 
tions have two periods rrii and m2- This is their essential property, (b) Our periodic and 
quasi-periodic models given in Sec. Ill are periodic or quasi-periodic in terms of the energy 
density p. (c) Furthermore, if the special values of the periods mi and m2 are infinite, these 
models in Sec. Ill are transformed into or have limits of the Chaplygin gas models, which 
can be seen from the formulae in (3.1). Thus, in this sense to reconstruct these models is 
considered to be two periodic or quasi periodical generalizations of the CG models. This is 
the justification of application of the Weierstrass functions to cosmological models, (d) Our 
models presented in Sec. IV are periodic on the cosmic time t (which is here dimensionless) 
contrast to the models in Sec. Ill which are periodic/quasiperiodic on p. (e) Moreover, our 
periodic/quasiperiodic models in Sec. Ill are singular at p = and those in Sec. IV at t = 0. 

It is also important to emphasize that the cosmological advantage to acquire the periodic 
evolution of the universe is to realize scenarios to avoid a Big Bang singularity, the finite-time 
future singularities and a Big Crunch singularity. By applying the obtained results and the 
discussed procedure to scenarios for the avoidance of singularities, e.g., the cyclic universe, 
the ekpyrotic scenario and the bouncing universe, we can find a non-singular cosmology. 

Furthermore, it has explicitly been illustrated that there are three type models with 
realizing the cosmological evolution of the EoS for dark energy : (i) The universe always 
stays in the non-phantom (quintessence) phase (the MG - XIV model), (ii) The universe 



20 



always evolves within the phantom phase (the MG - XII, MG - XIII, MG - XV, MG - XVI 
and MG - XVII models), (iii) The crossing of the phantom divide can be realized (the MG 
- XVIII, MG - XIX, MG - XX and MG - XXXIII models). If the universe always stays in 
the phantom phase, it is impossible to describe the whole evolution history of the universe, 
i.e., the decelerated phases such as the radiation and matter dominated stages. Thus, these 
models are ruled out. In addition, the corresponding description of a canonical scalar field 
model to the models in the ab ove (iii) category is not physical but just mathematical because 



in the light of quintom model [124l |. it is impossible for a single canonical scalar field to realize 



the crossing of the phantom divide. 

The scalar fields and the corresponding potentials have been analyzed for different types 
of MG models mentioned above. The EoS parameters have been derived and their natures 
have also been illustrated graphically during the evolution of the universe. 

On the other hand, at the current stage the cosmological constant is consistent with the 
observational data jj], but there still exists the possibility of dynamical dark energy, which 
is, e.g., described by a scalar field or a fluid. Thus, when we acquire the resul ts o f data 



analysis of the more precise future experiments like by the PLANCK satellite |l36l |. it is 
strongly expected that the investigations on the phase of the universe, i.e., the non-phantom 
or phantom phase, or the crossing of the phantom divide can become more meaningful. In 
the future work, it would be interesting that whether the above models may be suitable 
candidates for dark energy is tested by coming observational data fittings. 

Acknow^ledgments 

K.B. would like to sincerely acknowledge very kind and warm hospitality at Eurasian 
National University very much, where this work has been completed. In addition, he is 
grateful to National Center for Theoretical Sciences and National Tsing Hua University 
very much for the very kind and warm hospitality, where the revision of this work was 
executed. 



[1] D. N. Spergel et al. [WMAP Collaboration], Astrophys. J. Suppl. 148, 175 (2003) 
[arXiv:astro-ph/0302209] . 

21 



D. N. Spergel et al. [WMAP Collaboration], Astrophys. J. Suppl. 170, 377 (2007) 
|arXiv:astro-ph/0603449| . 

[3] E. Komatsu et al. [WMAP Collaboration], Astrophys. J. Suppl. 180, 330 (2009) 
far Xiv:0803.0"547] [astro-ph]]. 

E. Komatsu et al. [WMAP Collaboration], Astrophys. J. Suppl. 192, 18 (2011) 
[arXiv: 1001. 45381 [astro-ph.CO]]. 

[5] S. Perlmutter et al. [SNCP Collaboration], Astrophys. J. 517, 565 (1999) 

|arXiv:astro-ph/9812133 ) . 
[6] A. G. Riess et al. [Supernova Search Team Collaboration], Astron. J. 116, 1009 (1998) 

|arXiv:astro-ph/980520l] . 
[7] M. Tegmark et al. [SDSS Collaboration], Phys. Rev. D 69, 103501 (2004) 

|arXiv:astro-ph/0310723| . 
[8] U. Seljak et al [SDSS Collaboration], Phys. Rev. D 71, 103515 (2005) 

|arXiv:astro-ph/0407372 | . 
[9] D. J. Eisenstein et al. [SDSS Collaboration], Astrophys. J. 633, 560 (2005) 

|arXiv:astro-ph/050lT7T] . 



[10] B. Jain and A. Taylor, Phys. Rev. Lett. 91, 141302 (2003) |[arXiv:astro-ph/0306 046 1 . 

[11] E. J. Copeland, M. Sami and S. Tsujikawa, Int. J. Mod. Phys. D 15, 1753 (2006) 

|arXiv:hep-th/06030"57] . 
[12] R. Durrer and R. Maartens, Gen. Rel. Grav. 40, 301 (2008) [arXiv: 071 1.0077) [astro-ph]]. 
[13] R. Durrer and R. Maartens. rarXiv:0811.4T32l [astro-ph]. 
[14] Y. F. Cai, E. N. Saridakis, M. R. Setare and J. Q. Xia, Phys. Rept. 493, 1 (2010) 

|arXiv:0909.2776l [hep-th]]. 
[15] S. Tsujikawa, larXiv: 1004. 1493) [astro-ph.CO]. 

[16] L. Amendola and S. Tsujikawa, Dark Energy (Cambridge University press, 2010). 
[17] M. Li, X. D. Li, S. Wang and Y. Wang, Commun. Theor. Phys. 56, 525 (2011) 

|arXiv:1103.5"870] [astro-ph.CO]]. 
[18] K. Bamba, S. Capozziello, S. Nojiri and S. D. Odintsov, Astrophys. Space Sci. 342, 155 

(2012) jar Xiv: 1205.34211 [gr-ac]]. 
[19] S. Nojiri and S. D. Odintsov, Phys. Rept. 505, 59 (2011) |arXiv:1011.0"544] [gr-ac]]. 
[20] S. Nojiri and S. D. Odintsov, eConf C0602061, 06 (2006) [Int. J. Geom. Meth. Mod. Phys. 

22 



4, 115 (2007)] [ arXiv:hep-th/0601213] . 
[21] T. P. Sotiriou and V. Faraoni, Rev. Mod. Phys. 82, 451 (2010) [arXiv: 0805. 1726) [gr-ac]]. 
[22] S. Capozziello and V. Faraoni, Beyond Einstein Gravity (Springer, 2010). 
[23] S. Capozziello and M. De Laurentis, Phys. Kept. 509, 167 (2011) [arXiv:1108.6266l [gr-ac]]. 
[24] A. De Felice and S. Tsujikawa, Living Rev. Rel. 13, 3 (2010) garXiv:1002.4928l [gr-oc]]. 
[25] T. Clifton, P. G. Ferreira, A. Padilla and C. Skordis, Phys. Rept. 513, 1 (2012) 

|arXiv:1106.2476] [astro-ph.CO]]. 
[26] S. Capozziello, M. De Laurentis and S. D. Odintsov, Eur. Phys. J. C 72, 2068 (2012) 

larXiv: 1206.4842] [gr-qc]]. 
[27] R. R. Caldwell, M. Kamionkowski and N. N. Weinberg, Phys. Rev. Lett. 91, 071301 (2003) 

|arXiv:astro-ph/0302506| . 
[28] B. Mclnnes, JHEP 0208 (2002) 029 [arXiv:hep-th/0112066| . 



[29] S. Nojiri and S. D. Odintsov, Phys. Lett. B 562, 147 (2003) |arXiv:hep-th/0303117| . 



[30] S. Nojiri and S. D. Odintsov, Phys. Lett. B 571, 1 (2003) [arXiv:hep-th/0306212| . 

[31] V. Faraoni, Int. J. Mod. Phys. D 11, 471 (2002) |arXiv:astro-ph/01 10067] . 

[32] P. F. Gonzalez-Diaz, Phys. Lett. B 586, 1 (2004) |arXiv:astro-ph/0312579l . 

[33] E. Ehzalde, S. Nojiri and S. D. Odintsov, Phys. Rev. D 70, 043539 (2004) 

|arXiv:hep-th/04050"34] . 
[34] P. Singh, M. Sami and N. Dadhich, Phys. Rev. D 68, 023522 (2003) |arXiv:hep-th/0305110] . 



[35] C. Csaki, N. Kaloper and J. Terning, Annals Phys. 317, 410 (2005) |arXiv:astro-ph/0 409596] . 
[36] P. X. Wu and H. W. Yu, Nucl. Phys. B 727, 355 (2005) |arXiv:astro-ph/0407424| . 
[37] S. Nesseris and L. Perivolaropoulos, Phys. Rev. D 70, 123529 (2004) 
|arXiv:astro-ph/0410309| . 



[38] M. Sami and A. Toporensky, Mod. Phys. Lett. A 19, 1509 (2004) |arXiv:gr-qc/0312009|. 



[39] H. Stefancic, Phys. Lett. B 586, 5 (2004) |arXiv:astro-ph/0310904| . 

[40] L. P. Chimento and R. Lazkoz, Phys. Rev. Lett. 91, 211301 (2003) [arXiv:gr-qc/0307111] . 

[41] J. G. Hao and X. Z. Li, Phys. Lett. B 606, 7 (2005) |arXiv:astro-ph/0404154| . 

[42] E. Ehzalde, S. Nojiri, S. D. Odintsov and P. Wang, Phys. Rev. D 71, 103504 (2005) 

|arXiv:hep-th/05020"82] . 
[43] M. P. Dabrowski and T. Stachowiak, Annals Phys. 321, 771 (2006) |arXiv:hep-th/0411199] . 
[44] F. S. N. Lobo, Phys. Rev. D 71, 084011 (2005) [arXiv:gr-qc/0502099| . 

23 



[45 
[46; 

[47; 
[48; 
[49; 

[50 

[51 

[52; 
[53; 
[54; 

[55 
[56 

[57; 
[58; 

[59; 

[60 
[61 
[62 
[63 

[64; 

[65; 

[66 



R. G. Cai, H. S. Zhang and A. Wang, Commun. Theor. Phys. 44, 948 (2005) 

[arXiv:hep-th/0505186] . 

I. Y. Aref'eva, A. S. Koshelev and S. Y. Vernov, Phys. Rev. D 72, 064017 (2005) 

|arXiv:astro-ph/0507067l . 



H. Q. Lu, Z. G. Huang and W. Fang, arXiv:hep-th/0504038 



W. Godlowski and M. Szydlowski, Phys. Lett. B 623, 10 (2005) |arXiv:astro-ph/0507322] . 
J. Sola and H. Stefancic, Phys. Lett. B 624, 147 (2005) |arXiv:astro-ph/0505133| . 

B. Guberina, R. Horvat and H. Nikohc, Phys. Rev. D 72, 125011 (2005) 
|arXiv:astro-ph/0507666] . 

Y. Shtanov and V. Sahni, Class. Quant. Grav. 19, LlOl (2002) |arXiv:gr-qc/0204040| . 

J. D. Barrow, Class. Quant. Grav. 21, L79 (2004) |arXiv:gr-qc/0403084] . 

S. Nojiri and S. D. Odintsov, Phys. Lett. B 595, 1 (2004) [arXiv:hep-th/0405078| . 

S. Nojhi and S. D. Odintsov, Phys. Rev. D 70, 103522 (2004) |arXiv:hep-th/ 0408170) . 

S. Nojiri and S. D. Odintsov, Phys. Rev. D 72, 023003 (2005) | arXiv:hep-th/0505215| . 

S. Cotsakis and I. Klaoudatou, J. Geom. Phys. 55, 306 (2005) jarXiv:gr-qc/ 0409022]. 

M. P. Dabrowski, Phys. Rev. D 71, 103505 (2005) [arXiv:gr-qc/0410033| . 

L. Fernandez-Jambrina and R. Lazkoz, Phys. Rev. D 70, 121503 (2004) 

|arXiv:gr-qc/0410124] . 

L. Fernandez-Jambrina and R. Lazkoz, Phys. Lett. B 670, 254 (2009) [ arXiv: 0805. 2284) [gr- 

qc]]. 

J. D. Barrow and C. G. Tsagas, Class. Quant. Grav. 22, 1563 (2005) |arXiv:gr-qc/0411045] . 

H. Stefancic, Phys. Rev. D 71, 084024 (2005) [ arXiv:astro-ph/0411630| . 

C. Cattoen and M. Visser, Class. Quant. Grav. 22, 4913 (2005) |arXiv:gr-qc/0508045] . 

P. Tretyakov, A. Toporensky, Y. Shtanov and V. Sahni, Class. Quant. Grav. 23, 3259 (2006) 

|arXiv:gr-qc/0510104] . 

A. Balcerzak and M. P. Dabrowski, Phys. Rev. D 73, 101301 (2006) |arXiv:hep-th/0604034] . 

M. Sami, P. Singh and S. Tsujikawa, Phys. Rev. D 74, 043514 (2006) |arXiv:gr-qc/0605113] . 

M. Bouhmadi-Lopez, P. F. Gonzalez-Diaz and P. Martin-Moruno, Phys. Lett. B 659, 1 



(2008) [arXiv:gr-qc/0612135i . 
[67] A. V. Yurov, A. V. Astashenok and P. F. Gonzalez-Diaz, Grav. Cosmol. 14, 205 (2008) 
|arXiv:0705.4T08] [astro-ph]]. 

24 



[68] T. Koivisto, Phys. Rev. D 77, 123513 (2008) |arXiv: 0803 .33991 [gr-ac]]. 

[69] J. D. Barrow and S. Z. W. Lip, Phys. Rev. D 80, 043518 (2009) |arXiv:0901. 16261 [gr-ac]]. 

[70] M. Bouhmadi-Lopez, Y. Tavakoli and P. V. Moniz, JCAP 1004, 016 (2010) [arXiv:0911.1428l 

[gr-qc]]. 
[71] S. Nojiri, S. D. Odintsov and S. Tsujikawa, Phys. Rev. D 71 (2005) 063004 

|arXiv:hep-th/05010"25] . 
[72] S. Nojiri and S. D. Odintsov, Phys. Rev. D 78, 046006 (2008) |arXiv:0804.3"5T9l [hep-th]]. 
[73] K. Bamba, S. Nojiri and S. D. Odintsov, JCAP 0810 (2008) 045 |arXiv:0807.2575l [hep-th]]. 
[74] K. Bamba, S. D. Odintsov, L. Sebastiani and S. Zerbini, Eur. Phys. J. C 67 (2010) 295 

|arXiv:0911.4"390l [hep-th]]. 
[75] P. J. Steinhardt and N. Turok, Phys. Rev. D 65, 126003 (2002) [arXiv:hep-th/0111098| . 
[76] P. J. Steinhardt and N. Turok, |arXiv:hep-th/01 11030] . 



[77] J. Khoury, B. A. Ovrut, P. J. Steinhardt and N. Turok, Phys. Rev. D 66, 046005 (2002) 

|arXiv:hep-th/01090"50] . 
[78] P. J. Steinhardt and N. Turok, Nucl. Phys. Proc. Suppl. 124, 38 (2003) 

|arXiv:astro-ph/0204479| . 
[79] P. J. Steinhardt and N. Turok, Science 296, 1436 (2002). 
[80] J. Khoury, P. J. Steinhardt and N. Turok, Phys. Rev. Lett. 92, 031302 (2004) 

|arXiv:hep-th/0307132] . 
[81] P. J. Steinhardt and N. Turok, Science 312, 1180 (2006) |astro-ph/ 0605173]. 
[82] K. Saaidi, H. Sheikhahmadi and A. H. Mohammadi, Astrophys. Space Sci. 338, 355 (2012) 

|arXiv:1201.0"275l [gr-qc]] . 
[83] S. Nojiri, S. D. Odintsov and D. Saez-Gomez, AIP Conf. Proc. 1458, 207 (2011) 

|arXiv:l 108.07671 [hep-th]]. 
[84] Y. F. Cai and E. N. Saridakis, J. Cosmol. 17, 7238 (2011) [ar2£iYin08.6052" [gr-qc]]. 
[85] V. Sahni and A. Toporensky, Phys. Rev. D 85, 123542 (2012) [arXiv: 1203.0395 [gr-qc]]. 
[86] D. Y. Chung, | arXiv:physic7/010506l | 
[87] J. Khoury, B. A. Ovrut, P. J. Steinhardt and N. Turok, Phys. Rev. D 64, 123522 (2001) 

|arXiv:hep-th/01032"39] . 



[88] J. Khoury, B. A. Ovrut, P. J. Steinhardt and N. Turok, arXiv:hep-th/0105212 



[89] R. Y. Donagi, J. Khoury, B. A. Ovrut, P. J. Steinhardt and N. Turok, JHEP 0111, 041 

25 



(2001) [arXiv:hep-th/0105199| . 
[90] D. N. Page, Class. Quant. Grav. 1, 417 (1984). 



[91] P. Peter and N. Pinto-Neto, Phys. Rev. D 65, 023513 (2001) [arXiv:gr-qc/0109038| . 
[92] P. Peter and N. Pinto-Neto, Phys. Rev. D 66, 063509 (2002) |arXiv:hep-th/0203013| . 



[93] Y. Shtanov and V. Sahni, Phys. Lett. B 557, 1 (2003) [ arXiv:gr-qc/0208047| . 

[94] T. Biswas, A. Mazumdar and W. Siegel, JCAP 0603, 009 (2006) |arXiv:hep-th/0508194| . 

[95] Y. F. Cai, T. Qiu, Y. S. Piao, M. Li and X. Zhang, JHEP 0710, 071 (2007) |arXiv:0704.1090l 

[gr-qc]]. 
[96] P. Creminelh and L. Senatore, JCAP 0711, 010 (2007) |arXiv:hep-th/0702165| . 
[97] Y. S. Piao, Phys. Lett. B 677, 1 (2009) |arXiv:0901.2644] [gr-ac]]. 
[98] J. Zhang, Z. G. Liu and Y. S. Piao, Phys. Rev. D 82, 123505 (2010) [arXiv: 1007.2498) [hep- 

th]]. 
[99] Y. S. Piao, Phys. Lett. B 691, 225 (2010) j arXiv:1001.0"63ri [hep-th]]. 
[100] Z. G. Liu and Y. S. Piao, Phys. Rev. D 86, 083510 (2012) jar Xiv: 1201. 1371] [gr-ac]]. 
[101] M. Novello and S. E. P. Bergliaffa, Phys. Rept. 463, 127 (2008) |arXiv:0802. 16341 [astro-ph]]. 
[102] R. Myrzakulov, to appear in Gen. Rel. Grav. DOI: 10.1007/sl0714-012-1439-z, 

larXiv:1008.4486l [physics.gen-ph]. 
[103] R. Myrzakulov, Advances in High Energy Physics, 2012, 868203 (2012) [arXiv:1204.1093l 

[physics.gen-ph]]. 
[104] K. Esmakhanova, Y. Myrzakulov, G. Nugmanova and R. Myrzakulov, Int. J. Theor. Phys. 

51, 1204 (2012) [arXiv: 1102.4456) [phvsics.gen-ph]]. 
[105] K. R. Yesniakhanova, N. A. Myrzakulov, K. K. Yerzhanov, G. N. Nugmanova, N. S. Serik- 

bayaev and R. Mvrzakulov. larXiv:1201. 4360 1 [phvsics.gen-ph]. 
[106] G. W. Gibbons and M. Vyska, Class. Quant. Grav. 29, 065016 (2012) [arXiv:1110.6508) 

[gr-qc]]. 
[107] L Bochicchio, S. Capozziello and E. Laserra, Int. J. Geom. Meth. Mod. Phys. 8, 1653 (2011) 

|arXiv:1111.08T9) [gr-qc]] . 
[108] B. G. Dimitrov, J. Math. Phys. 44, 2542 (2003) |arXiv:hep-th/0107231| . 
[109] J. D'Ambroise,|ar2£iYiQ908J481' [gr-qc]. 

[110] M. Bouhmadi-Lopez, L. J. Garay and P. F. Gonzalez-Diaz, Phys. Rev. D 66, 083504 (2002) 
p[?XiV:gr-qc/0204072] . 

26 



[Ill] K. Bamba, K. Yesmakhanova, K. Yerzhanov and R. Mvrzakulov. larXiv:1203.3401l [gr-qc], 
[112] A. Y. Kamenshchik, U. Moschella and V. Pasquier, Phys. Lett. B 511 (2001) 265 

[ailGvigr-qc/OlOSOOi] . 
[113] M. C. Bento, O. Bertolami and A. A. Sen, Phys. Rev. D 66 (2002) 043507 

|arXiv:gr-qc/0202064] . 



[114] H. B. Benaoum, arXiv:hep-th/0205140 



[115] D. Saez-Gomez, Grav. Cosmol. 15, 134 (2009) ^ arXiv:0804.4586] [hep-th]]. 

[116] M. Sharif, K. Yesmakhanova, S. Rani and R. Myrzakulov, Eur. Phys. J. C 72, 2067 (2012) 

|arXiv:1204.2T8T] [physics.gen-ph]]. 
[117] K. Bamba, S. Nojiri, S. D. Odintsov and M. Sasaki, Gen. Rel. Grav. 44, 1321 (2012) 

|arXiv:1104.2692] [hep-th]]. 
[118] U. Alam, V. Sahni and A. A. Starobinsky, JCAP 0406, 008 (2004) |arXiv:astro-ph/0403687l . 
[119] U. Alam, V. Sahni and A. A. Starobinsky, JCAP 0702, Oil (2007) |arXiv:astro-ph/061238l| . 



[120] S. Nesseris and L. Perivolaropoulos, JCAP 0701, 018 (2007) [arXiv:astro-ph/0610092] . 

[121] P. U. Wu and H. W. Yu, Phys. Lett. B 643, 315 (2006) [arXiv: astro-ph/06 1 1507| . 

[122] H. K. Jassal, J. S. Bagla and T. Padmanabhan, Mon. Not. Roy. Astron. Soc. 405, 2639 

(2010) |larXiv:astro-ph/0601389] . 
[123] Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables 

(1964), edited by M. Abramowitz and 1. A. Stegun. 
[124] G. -B. Zhao, J. -Q. Xia, M. Li, B. Feng and X. Zhang, Phys. Rev. D 72, 123515 (2005) 

[astro-ph/0507482| . 
[125] T. Chiba, T. Okabe and M. Yamaguchi, Phys. Rev. D 62, 023511 (2000) 



|arXiv:astro-ph/9912463 1 . 



[126] C. Armendariz-Picon, T. Damour and V. F. Mukhanov, Phys. Lett. B 458, 209 (1999) 

[aFxiv:hep-th/9904075] . 
[127] J. Garriga and V. F. Mukhanov, Phys. Lett. B 458, 219 (1999) |arXiv:hep-th/9904176] . 
[128] C. Armendariz-Picon, V. F. Mukhanov and P. J. Steinhardt, Phys. Rev. Lett. 85, 4438 (2000) 

[M^v:astro-ph/0004134 | . 
[129] C. Armendariz-Picon, V. F. Mukhanov and P. J. Steinhardt, Phys. Rev. D 63, 103510 (2001) 



|arXiv:astro-ph/0006373] . 



[130] R. de Putter and E. V. Linder, Astropart. Phys. 28, 263 (2007) [arXiv:0705.0400] [astro-ph]]. 

27 



[131] A. Nicolis, R. Rattazzi and E. Trincherini, Phys. Rev. D 79, 064036 (2009) |arXiv:081 1.21971 

[hep-th]]. 
[132] C. DefFayet, G. Esposito-Farese and A. Vikman, Phys. Rev. D 79, 084003 (2009) 

fa rXiv:0901.13T4] [hep-th]]. 
[133] C. Deffayet, S. Deser and G. Esposito-Farese, Phys. Rev. D 80, 064015 (2009) 

|arXiv:0906. 19671 [gr-qc]] . 
[134] C. Deffayet, S. Deser and G. Esposito-Farese, Phys. Rev. D 82, 061501 (2010) 

|arXiv:1007.5278 [gr-qc]]. 
[135] N. Shirai, K. Bamba, S. Kumekawa, J. Matsumoto and S. Nojiri, Phys. Rev. D 86, 043006 



(2012) jar Xiv: 1203.49621 [heo-th]]. 



[136] http://www.sciops.esa.int/index.php?project=PLANCK 



28 



